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Abstract 

By using a parametrization of the non-linear Walecka model which takes into account the binding 
energy of different hyperons, we present a study of particle production yields measured in central 
Au-Au collision at RHIC. Two sets of different hyperon-meson coupling constants are employed 
in obtaining the hadron production and chemical freeze-out parameters. These quantities show a 
weak dependence on the used hyperon-meson couplings. Results are in good overall accordance 
with experimental data. We have found that the repulsion among the baryons is quite small and, 
through a preliminary analysis of the effective mesonic masses, we suggest a way to improve the 
fittings. 

PACS numbers: 21.65.-f,24.10.Jv,25.75.-q 
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I. INTRODUCTION 

Lattice simulations of QCD indicate that, at zero baryon number density and temper- 
atures of the order of 150-170 MeV, quarks and gluons become deconfined [l|, y, yj. This 
novel phase of nuclear matter is commonly called quark-gluon plasma (QGP). It is believed 
that in nature temperatures of the order of 150 MeV existed only shortly after the Big Bang. 
For more than two decades, attempts have been made to recreate similar conditions within 
the collision of heavy nuclei at ultra-relativistic energies, looking for signatures of the pro- 
duction of the QGP, which subsequently hadronizes. The hadronization stage has been well 
discussed in the literature j^. After this process takes place, one assumes that a hadronic 
gas appears in a sizeable region of space. Hydrodynamic models were used to describe this 
hadronic evolution, based on the hypothesis that the fireball size is large in comparison 
with the mean free path of the hadrons [5|, l6|. One may ask whether the hydrodynamic 
approach needs to consider viscosity and heat conduction effects, even before the freeze-out 
occurs. Here, one must say that the employment of any kind of model has to be taken with 
caution. In the gas regime, and for temperatures far from the system critical temperature, 
the relaxation times which control the chemical equilibrium of hadrons may become large 
compared with collision times, causing the gas to reach the chemical equilibrium at an early 







stage of the expansion, favoring an overpopulation of pions [7|. Temperature and chemical 
potential are related to the chemical freeze-out regime. At low beam energies the chemi- 
cal freeze-out of the resulting gas may be described by low temperature and high baryonic 
chemical potential. As the beam energy increases this situation reverts. The mesons become 
more numerous, the temperature increases and the chemical baryonic potential decreases. 
For small baryonic chemical potential and high temperature the number of baryons becomes 
close to the anti-baryons. Within this scenario, a QGP-hadronic matter phase transition is 
probably going to happen in heavy-ion collisions at very high energy, as the ones extected 
to take place soon at the Large Hadron Collider (LHC/ ALICE). 

Hadron multiplicities are observables which can provide valuable information on the na- 
ture of the medium from which they are produced. In this way, they are natural tools to 
look for QGP signatures. Different models have been used to improve the understanding 
of the QGP-hadronic matter phase transition [8|, |9|, [lO|, lul, 112] ■ The use of thermal models 
to describe hadronic collision spectra has also been considered |l3|, Uj^ . Curiously, despite 



the fact that no microscopic theoretical basis was given for thermahzation in such colhsions, 
these models have been quite successful. 

In thermal models, the statistical distribution is controlled basically by the temperature 
and the chemical potential of a free hadron-gas. However, these ingredients alone does not 
suffice since without a repulsion dynamics, the gas does not expand properly. Therefore, 
for baryons and mesons, an excluded volume is needed in the same spirit as in a Van der 
Waals gas. In particular, a detailed explanation of how incorporate a thermodynamically 
consistent excluded volume effect in a hadron gas was presented in 
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16[ |. The excluded 



volume parameters are determined accordin g to the fits of heavy-ion collision data. To fit 



13l . Il4l . [iTJ], a common hard-core for baryons and 



Au-Au and Pb-Pb hadronic ratios data 
mesons suffices. 

In a different approach, that does not consider any correlation between baryons, an 
interesting work 18j has shown that the energy per baryon remains almost constant and 
equal to 1 GeV in a broad range of chemical freeze-out temperatures. 

More recently, studies including hadron-hadron dynamics in a sophisticated relativistic 
chiral SU(3) model, have shown that the hadron ratios obtained from SIS, AGS, SPS to 

n 

RHIC hadron energies are also well described [19|. This calculation, contrary to thermal 
models, takes into account in medium hadronic mass effects. This new ingredient appears to 
modify substantially the temperature and the baryonic chemical potential to the freeze-out 
fittings. This happens even when the SU(3) models are undergoing a first or a second order 
phase transition. As a side information, let us remark that a first order phase transition in 
a hadronic model induces a dramatic decrease in the effective baryonic masses as well as a 
discontinuity in the entropy. In the SU(3) model, both mesonic and baryonic masses are 
obtained with in medium effects. 

Between the simplicity of the thermal models 
ativistic chira 
models 
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17( 1 and the complexity of the rel- 
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SU(3) model [l9|], one would ask whether well known relativistic hadronic 
211] are able to successfully fit the mentioned chemical freeze-out hadronic ratios. 
In a previous study [2^] the antibaryon (p. A) production in relativistic nuclear collisions us- 
ing nonlinear models has been presented. The authors estimated the ratio of antibaryon and 
baryon densities in interacting and in free systems. It was also pointed out that the anni- 
hilation of antibaryons in surrounding matter at the final stage of the reaction may reduce 



their abundancy 221]. Relativistic hadronic models are based on field theories describing 



nuclear matter as a strongly interacting system of baryons and mesons. This approach is 
assumed to be valid below the deconfinement phase transition, and is based on the identi- 
fication of the appropriate degrees of freedom at this scale. The prototype of such theory 



is Quantum Hadrodynamics (QHD) 20|, |2l|, which has been shown to describe wel. 
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25,^ 



many 



27|. This 



properties of nuclear matter, finite nuclei and neutron stars |20l . 
theory can be extended to include many-body correlations as density-dependent meson cou- 
plings [28|, l29|] . The successful description of nuclear properties indicates that the essential 
aspects of low energy strong interactions are well described by QHD. At high energy, it 
can describe dynamically the evolution from the hadronic fireball starting scenario to the 
chemical freeze-out in which the baryons and the light mesons populate the gas. One of 
the advantages of using relativistic quantum-hadronic models is that they can also describe 



the already mentioned QGP-hadronic matter phase transition as the chemical potential in- 



creases 
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12l | . The physics of the hadronic models at high temperature is very rich and 



its behavior has been well investigated for zero chemical potential 30| . Depending on their 
parametrization, the models can present a first or second order phase transition in the same 
way as the already mentioned SU(3) model for the hadronic masses. 

It is natural then to ask whether these routinely used relativistic hadronic models can 
describe the hadron production observed today in Au-Au collisions at RHIC [ij]. Recently, 
some of these models, with and without density depend hadron-meson coupling constants, 
have been used in order to answer that question 311] • In that calculation not only the bary- 
onic octet was included but also the decuplet. A good overall description of the particle 
yields data was obtained and in particular for K^* /h~ and K^*/h~ ratios, a good improve- 
ment was achieved in comparison with the results obtained in the free gas approximation 
14| . In 31| it was shown that the new dynamics achieved by the models with density depend 
parameters makes the difference by having the best fits and that the mesons play an impor- 
tant role. To go further with this investigation, we have chosen to consider in the present 
work only the baryonic octet because the inclusion of the decuplet produced only a small 



change in the overa 



of the hyperons 32, 



1 results. Moreover, we take into account the different binding energies 
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35| in obtaining the chemical freeze-out parameters. We also 
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24|. 



restrict ourselves to one parametrization of the non-linear Walecka model [20|, 

In order to investigate the role played by the mesons, we have introduced, in a very crude 
and naive way, an ad-hoc effective mesonic mass m* that increases in the medium and obeys 



a universal scaling with its bare value: m* /m = a, a > 1, for all mesons. Using a best fitting 
analysis, we established a correlation between the x^ psr degree of freedom {x\of) ^^^ a, 
understanding the result as an indication that the increase of the in-medium meson masses 
introduces additional repulsion among the particles and mimics the repulsive effect obtained 
in other approaches, where the hadron excluded volume is taken into account 13|, [ij] . 

The paper is organized as follows. In Sec. II we present the models with a detailed 
discussion on chemical, strangeness and charge equilibrium. In Sec. Ill, we present the 
parametrizations used for the couplings between hyperons and mesons. In Sec. IV, we 
present our results and summarize our main conclusions. 

II. THE MODEL 

We assume that the chemical freeze-out can be described as a mixture of the lightest 
baryons and mesons. We use the following Lagrangian 

^ = ^QHD + ^irKpi (1) 



where 



with 



^QHD 



^b + ^m, (2) 



j 
1 



+ -W?„UJ^U^ - 



-Ad,pi-d^plY + \mipy^. 



Mesons {-n^ ,K^ ,K^ ,K^ , p^ ,K*^ ,K*'^ ,K*^^ enter Lagrangian Lt,kp in ([T]) as a free gas of 
bosons. Lagrangian ([2]) has a general form commonly used in relativistic mean field models. 
The sum over j in ([3]) extends over the octet of lightest baryons {n, p. A, S", S°, S+, S~, H"}, 
and their antiparticles. /3J is the corresponding isospin quantum number. The baryon and 



meson masses are listed in Table [B The values for the couplings are taken from 26| and 
shown in Table [ITl The Euler-Lagrange equations of motion for the mediating mesons in 
a system with translational and rotational invariance, within the mean field approximation, 



A^A S E TT K p K* 
m (MeV) 939 1116 1193 1318 138 495 776 893 



TABLE I: Masses used in the calculation. 



K (MeV) m*/m l„ (fm-^) /^ (fm-^) l^ (fm-^) b {im.-^) c 

300 0.70 11.79 7.149 4.411 0.01402 -0.001070 

TABLE II: Parameters used in the calculation (/j = gfj^/mf), taken from Ref. [26] {b here includes 
the product with mjv in Eq. (1) of this reference). 



are 

K = / J E ^'^J^sj - hV^ - cV^ J , (4) 

j 

Vp = lpY.^pJ^3J^ (6) 

j 

with K = gaN(^, K; = gujNUJQ, Vp = QpNPo, Xij = Qij / QiN and k = gf^/^h i = cr,^,P- The 
densities entering the above equations are 

1 r m* 
n 



t Tfl ■ 



^J 



1 

vr- 
n^j = hjUj, (9) 



^i = ^ / (fi+ - fj-)P^dp, (8) 



where m* = rrij — g^^ja is the effective mass of the baryon j, and ej{p) = Jp^ + m*, . The 
baryonic distribution functions fj± are the Fermi-Dirac ones 

^'^ = exp[(e,T^,)/T] + l' ^^°^ 

where the sign minus (plus) accounts for particles (antiparticles). The baryonic effective 
chemical potential Vj is defined as 

^3 = Pj - 9i^j^o - gpjhjPo^ (11) 



being fij the therinodynainical chemical potential of baryon j. The particle density of the 
free mesons it, K, p and K* are calculated using the Bose- Einstein distribution function in 
the form 

"'=2^/ exp[(e,-!.0/T]-/ ^^' / = vr, iT, p, TT* (12) 



where ei = Jp^ + mf , pi is the chemical potential, and 7/ is the corresponding spin degen- 
eracy. 

The problem is solved in the grand canonical ensemble, imposing the constraints of baryon 
number, strangeness and electric charge conservation, given respectively by [13|: 

VY^n.Bi = Qb = Z + N, (13) 

i 

VY^n^S, = Q5 = 0, (14) 

i 

Vj2n,h, = Qj, = {Z-N)/2, (15) 

i 

where i runs over all particles in the system, and {Bi, St, I^i} are the baryonic, strangeness 
and isospin quantum numbers of particle i. Each conserved charge has a conjugated chemical 
potential, namely {QB,Qs,Qh} ^^ {f^B, f^s, f^h}- The temperature T and the baryochem- 
ical potential fiB are the two independent parameters of the model, while the volume of 
the fireball V, the strangeness chemical potential /i^, and the isospin chemical potential /i/3 
are fixed by the three constraints (TT3 l) -( TT5l) . The chemical potential of each particle in the 
system can be written as a linear combination of these three conjugate chemical potentials 
in the form 

fii = BifiB + hifJ'h + SifJ'S- (16) 

In the baryonic sector, a zero chemical potential implies a zero baryonic density since 
the number of baryons is equal to the anti-baryons. For the vr'', for instance, zero chemical 
potential is obtained directly from Eq. flT6l) . using i?^o = Ja^o = 5*7^0 = 0. Despite the fact 
that Pt^o = 0, the density for vr^ is different from zero. This remark applies only to vr'' and 
p°. In particular, 7r° results completely decoupled from the other particles. Its population 
is that of a boson gas at temperature T and zero chemical potential. This is the reason 
why it was not considered in Lagrangian Cj^xp in (III)- Consistently, particle ratios do not 
involve the vr" and p° mesons [14]. Hence, these particles do not need to be considered in 
the equations of the chemical potentials. Note that the p° meson appears in Lagrangian Cm, 
Eq. ([3]), only as a mediating field in the mean field approximation. 
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In the bosonic sector, the chemical potentials have to obey the physical inequality, 

IJ^i<fni, (17) 

where mi is the bosonic mass. This constraint is imposed naturally by statistical mechanics 
in order to keep the bosonic partition function finite. 

The constraints contained in Eqs. (11611171) allow us to restrict the possible values of the 
generalized chemical potentials by the following inequalities: 

0<\fii,\<m^, (18) 

< l/isi < -^\l^h\ + rriK- (19) 

This set of inequalities was implemented numerically in order to solve the large self- 
consistent set of equations of the problem. 

III. THE HYPERONIC COUPLINGS 

The coupling constants of the hyperons to the uj and p mesons are fixed using SU(6) 



symmetry [35|, which means that the vector coupling constants scale with the number of 
light quarks in the baryon in the form g^^ : g^^ : g^^, : g^^ = 3:2:2:1. The ratios of the 
vector couplings are then fixed to 

_ 2 _ 1 

meaning that the vector coupling constants of the hyperonic sector are fixed once the gt^N 
coupling constant is know. 

The coupling strength of the p meson, however, is given by the isospin of the baryon in 
the form gp^ : gp\ : gps : gps = 1/2 : : 1 : 1/2. Again, the isovector coupling constants of 
the hyperonic sector are fixed once gpN is know. This is taken into account automatically 
with the specific form of the isovector couplings in Lagrangian ([3]) choosing 

Xpi = l, ^ = A,S,H. (21) 

The couphng constants {xa-j}i=A,E,H of the hyperons with the scalar meson a are constrained 
jy^the hyp ernuclear potentials in nuclear matter to be consistent with hypernuclear data 



32 



33 



3J]. The hypernuclear potentials were constructed as 

Vi = x^iV^N - XaiVaN, i = A, S, H, (22) 



•^ujA. •^cjS -^cjH -^(tA '^aYi ^a^ 



^pi 



Setl 0.6666 0.6666 0.6666 0.6104 0.6104 0.6104 0.6104 
Set2 0.6666 0.6666 0.3333 0.6106 0.4046 0.3195 1.000 



m- 



TABLE III: Parametrizations used for the hyperon coupling constants. Setl is taken from Ref. 
Set2 is obtained from Setl imposing SU(6) symmetry and different binding energies for hyperons 
in nuclear matter (see text for details). The last column, Xpj, refers to the coupling of the p meson 
with all hyperons, z = A, S, H. 



where V^n = 215.83 MeV and V^]^ = 281.47 MeV are the nuclear potentials for 



nuclear matter at saturation with the parameters of Table [TTl Following Ref. 33| we use 



symmetric 



Va = -28 MeV, V^ = 30 MeV, I^h = -18 MeV. (23) 

All hyperon coupling ratios {x^i, x^i, x pi}i=\^-£^s are now known once the coupling constants 
{QaN^QuN^QpN} of the nucleon sector are given. 

We have to guarantee also that the hyperon coupling ratios are constrained to the known 



neutron star maximum masses, as in Ref. 261]. In this reference windows are established 



for the values of x^ and x^^: x^ and x^i are the couplings between hyperons and a and u) 
meson respectively and they have to be compatible with the allowed maximum neutron star 
masses and with the A binding energy in nuclear matter. Couplings x^ and x^; are the same 
for all hyperons and the Xp couplings are fixed arbitrarily as Xp = x^j. The SU(6) value 
Xt^ = 2/3 of Eq.( l20l) together with x^ = 0.6104, is one of the allowed couple of values (they 



were interpolated from Table I of Re: 
mil The other couplings used in Ref. 



261]). This choice corresponds to Setl row in Table 

261 ] are those shown in Table [TTl Through this work we 

refer to this relativistic model as GMl. Now we modify Setl to take into account constraints 

( I20|) - (l23|) . The value x^a = 2/3 is kept constant and the others are modified accordingly 

with the constraints. The resulting set of hyperonic couplings is the Set2 row in Table IIIII 

To see the impact of this redefinition of hyperonic couplings on neutron star properties. 



we show in Fig. [T]the equations of state calculated using the model of Ref. [26|] with Setl 
(dashed line) and with Set2 (full line) of Table [TTTl for the hyperon coupling constants. As we 
can see, both curves are almost indistinguishable, ensuring the same neutron star properties. 
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2x10 3x10 

£ [MeV/fm"^] 



FIG. 1: Equations of state for the relativistic models of Ref. [26i] (dashed line, Setl in Table 
imp and the one of this work (full line, Set2 in Table IIIip . Both equations of state are almost 
indistinguishable, ensuring the same neutron-star maximun masses. 

This validates Set2 as a set of hyperonic couplings that preserves maxiinuin neutron star 
masses as well as experimental results for the binding energy of hyperons in symmetric 
nuclear matter. In the following, we use Setl and Set2 for the hyperonic couplings in the 
Lagrangian ([1]) together with parameters of Table [TT] to investigate particle production in 
relativistic heavy ion collisions. 



IV. RESULTS AND CONCLUSIONS 

We apply the above formalism to the description of the experimental data for hadron pro- 
duction yield in Au-Au collisions at RHIC from STAR, PHENIX, PHOBOS and BRAHMS 
collaborations, following [iJ] and references therein. We adjust the free parameters T and 
fiB to get the best description of the data, based on a x^ analysis of the form [l3| 

/ pex'p. pmod. \ 2 

- ^ R, -R, ^ _ ^^^j 



X 



CTi 



where K^""^' is the ith particle ratio obtained from the models used here and Uj is the 
experimental error. When more than one experimental data is available, we calculated a 
weighted average. This leaves us with seven effective degrees of freedom for the fit. 

Our results are summarized in Table HVl for the hyperon coupling constants given in Table 
mil We also show, in the fourth column, the results from a non-interacting gas of baryons 

11 



GMl Free Free* Th Exp. Data Exp. 
Ratio Setl Set2 

p/p 0.673 0.676 0.674 0.638 0.629 0.65±0.07 STAR 

0.64±0.07 PHENIX 

0.60±0.07 PHOBOS 

0.64±0.07 BRAHMS 

p/tt- 0.048 0.046 0.038 0.080 0.078 0.08±0.01 STAR 

7r-/7r+ 1.005 1.005 1.004 1.011 1.007 1.00±0.02 PHOBOS 

0.95±0.06 BRAHMS 

K-/K+ 0.957 0.959 0.964 0.888 0.894 0.88±0.05 STAR 

0.78±0.13 PHENIX 

0.91±0.09 PHOBOS 

0.89±0.07 BRAHMS 

K-/-K- 0.239 0.237 0.231 0.175 0.145 0.149±0.02 STAR 

K*^/h- 0.063 0.062 0.059 0.038 0.037 0.06±0.017 STAR 

K*^/h- 0.060 0.059 0.057 0.033 0.032 0.058±0.017 STAR 

A/A 0.692 0.694 0.699 0.717 0.753 0.77±0.07 STAR 

S+/H- 0.723 0.713 0.725 0.806 0.894 0.82±0.08 STAR 
T (MeV) 148.6 148.0 145.7 169.4 174 
fiB (MeV) 33.27 32.51 28.88 38.52 46 
xlof 5.65 5.77 6.24 1.08 0.81 

TABLE IV: Comparisons of experimental particle ratios with the ones obtained from the relativistic 
mean-field models used in this work for the hyperon couplings belonging to Setl and Set2 , together 
with the chemical freeze-out temperature T, baryonic potential hb and x^^t- Also shown are the 
results from the free gas (Free), from the fr ee g as with ad-hoc effective meson masses (Free*, 



m* /m = 1.3), and from the thermal model of 



14l | (Th). Experimental values of Au-Au collisions 



at y/s = 130 MeV were taken from that reference and references therein. 
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and mesons (Free), in the fifth column the results from a non-interacting gas of baryons and 
mesons with effective meson masses as discussed below (Free*), together with the results 
from the thermal model of Ref. [ij] in column six (Th). 

We can see that both set of hyperonic couplings give almost the same results for the 
particle ratios. They are not very different from the model with zero coupling constants 
designated by Free, although in this case a smaller temperature and chemical potential 
were obtained. As a whole the fractions are well described, but we observe that there is a 
systematic deviation for the fractions involving mesons except for 7r^/7r+. This reflects the 
naive way the mesons were included. However, we point out that the K* //i~ and K* /h~ 



ratios have im pro ved with respect to the thermal model resu 



previous work 



ts in accordance with our 
For h~ we designated 



3l| and also with the chiral SU(3) calculations 
all negatively charged hadrons produced, including mesons and baryons. 

Regarding the quality of the fit, the resulting Xdof ^^ both sets are around 5.70, showing 
a slight improvement of Setl over Set2 . When these values are multiplied by the seven 
degrees of freedom used in our fitting, x^ reaches 40, which is a high value. This does not 
happen with the fits obtained by thermal models [13|, um, [17|] and with the chiral SU(3) 
model [l9(], for which x^ is around 6. 

Looking now at the freeze-out temperatures, we see that they are almost the same for 
the two sets of hyperonic couplings considered (T = 148.6 — 148.0 MeV), showing a small 
decrease when going from Setl to Set2 . This situation repeats for the chemical potential 
fitting (/i = 33.27 — 32.51 MeV). It is important to mention at this point that in these models 
the hadronic phase transition regions occur at temperatures T > 180 MeV [12], higher than 
the freeze-out temperature we found here. 

In order to better understand why we obtain a high Xdof^ ^^ argue that the repulsion 
among the baryons is quite small and cannot describe so well the chemical freeze-out envi- 
ronment. This happens because the chemical potential is too low in comparison with the 
temperature. The other cited approaches certainly have sufficient built in repulsion to avoid 
an overpopulation of particles. In the thermal model, it is done through excluded volumes 
introduced in a thermodynamically consistent way. On the other hand, the chiral SU(3) 
model includes meson-meson and meson-baryon interactions, providing effective masses for 
all hadrons, giving a higher kinetic energy which mimics an additional repulsion among the 
particles. 

13 



We think that the relativistic models themselves suffer from a lack of repulsive correlations 
at low densities, where the sources for the mediating meson fields are very weak. For example, 
in the GMl model with Set2 , studied in this work, we have Vo-n = 21.44 MeV, V^^ = 3.359 
MeV, and VpN = —0.8919 x 10^^ MeV. This fact results in a system almost identical to the 
one described by a free gas of baryons and mesons, as we can see comparing columns 2, 3 and 
4 in Table IIVI Therefore, the hadron production is not very sensitive to the hadron-meson 
interactions, manifested in the almost model independence of the Xdof obtained and also in 
the small dependence of the freeze-out parameters. In spite of the reduced number of models 



studied here, the results obtained in |3l|, where more models were considered, also reinforce 
this point. This conjecture may be extended to all relativistic hadronic models commonly 
used in the literature. 

Understanding the reason underlying the lack of repulsion in the description of the chem- 
ical freeze-out regime is a relevant question. In the following, we present a crude argument 
for that. Let us suppose, for simplicity, that baryon-baryon interaction proceeds through 
an attractive scalar {S) and a repulsive vector {V) meson-exchange potentials. In the mean 
field approach, at zero temperature and at normal saturated nuclear matter, the depth of 
S = 5 + V" is around 60 MeV. The temperature increase in these models favors the ap- 
pearance of anti-baryons. Table HVl shows that the p/p ratio is around 0.67. However, from 
G-parity symmetry, the anti-baryons themselves revert the sign of the vector potential in 



37 



such a way that H = S — V 36 



38|. Thus, anti-baryon-baryon interactions themselves 
do not carry repulsion. Therefore, the baryon-baryon repulsion alone has to compensate for 
the large attraction among the anti-baryons, explaining the lack of repulsion in the hadronic 
models. 

One way of implementing a repulsive content in this class of hadronic models would be 
by introducing an excluded volume for the particles in the system 15|, ll6|, l39(]. Of course, 
this approach requires the recalculation of the coupling constants in order to maintain the 
correct normal nuclear matter properties of the original models. By doing so, however, the 
hadronic models would loose its strong appealing due to their theoretical basis together with 
the success in reproducing so many observables in several different regimes of nuclear matter. 
Another possibility, inspired by the SU(3) hadronic model [l9|], is to include meson-meson 
self-coupling and meson-baryon couplings for the heavy mesonic sector, as K and p for 



example. We have studied this effect in an ad-hoc and crude way. In the statistical thermal 
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1.3 1.4 

m*/m 



FIG. 2: Xdof ^^ ^ function oim* /m (the same for all mesons), where m is the meson bare mass, for a 
free relativistic gas of baryons and mesons (baryon masses are mantained fixed to the experimental 
values of Table H). 

distributions of the free gas, we have artificially assumed in-medium effective mesonic masses 
m*, instead of the bare ones in vacuum. In this preliminary study, we have used the same 
value of m* /m for all mesons. In Figure [2] we see how x^^r becomes smaller as ■m*/m 
increases from 1, reaches a minimum around m* /m ^ 1.3, and starts to grow again. This 
shows that an increase in the effective in-medium mesonic mass goes in the right direction 
towards obtaining a better fit of the data. This is confirmed looking at column 5 in Table 
llVt where the particle ratios for m*/m = 1.3, close to the minimum of Fig. [2], are shown. 
The increase of the in-medium meson masses would be compatible with the notion of a 
thermal mass which increases linearly with temperature. This kind of behaviour can be seen 



in chiral effective models 



191] and also in the Walecka model 401]. Including the effective 



meson masses in this ad-hoc way, we automatically improve the ratios involving mesons, and 
reduce the fractions involving h~ . We point out, however, that in fact the uncertainty on the 
measured values of K*^ /h^ and K* /h~ are quite large. Also, the equilibrium temperatures 
increases more than 15% and the chemical potential 33%, coming closer to the results of the 
thermal model. 

This finding strongly suggests that introducing an appropriate dynamic treatment of the 
mesons, allowing for the modification of their masses in medium, would favor the correct 
data ratios without the need of including excluded volumes. However, the model dependence 



addressed in Ref. 



19| remains to be further investigated. 
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Another remark regards the particle ratios in the medium as we have calculated and the 
ones measured. The later comes from asymptotic particle states with vaccuum masses, i.e., 
on the mass shell {m* = m). In our model calculation the particles are off-shell and any 
mechanism to bring them on-shell would modify the ratios. However, we did not include 
annihilation and nonequilibrium processes which are claimed to reduce the calculated yields 



22|. 



Summarizing, in this work we have investigated the effects of different hyperon-meson 
coupling constants in hadron production and chemical freeze-out parameters. We have ob- 
served that the freeze-out temperatures and chemical potentials show a small dependence on 
the different parametrizations. Our results indicate that the repulsion among the baryons 
is quite small and the hadron production is not very sensitive to the hadron-meson inter- 
actions. A simple ad hoc in-medium mesonic mass study shows the importance of having 
meson-meson interaction in the hadronic models. 
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